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ABSTRACT 

We present the first high-resolution near-infrared images of the edge-on silhouette 
circumstellar disk, Orion 114-426, made using NICMOS on the Hubble Space Telescope. 
Images taken against the bright nebular background of the ionized hydrogen Paa line at 
1.87/xm show the major axis of the disk to be approximately 20% smaller than at 0.6/im, 
from which we deduce the structure of the edge of the disk. Continuum images of diffuse polar 
lobes above and below the plane of the disk show a morphology and evolution with wavelength 
consistent with predictions for reflection nebulae in a diffuse envelope with large polar cavities, 
surrounding a thin, massless, Keplerian disk, centered on an otherwise hidden central star. 
We make use of our observations and reasonable assumptions about the underlying disk 
structure to show that the disk mass is at least 10 M® and plausibly > 5 x 10~ 4 M Q . 

Subject headings: accretion, accretion disks — circumstellar matter — ISM: individual (Orion 
Nebula) — stars: formation, pre-main sequence — infrared: ISM: continuum, lines and bands 

1. Introduction 

The discovery of a family of circumstellar disks seen as dark silhouettes in projection against the 
Orion Nebula using the Hubble Space Telescope provided strong confirmation of the disk paradigm of 
star formation (O'Dell etal. 1993; O'Dell & Wen 1994; McCaughrean & O'Dell 1996 [M096]). The 
disks range in diameter from 50-1000 AU, and thus the ~50AU (0.1 arcsec at 450 pc) resolution of the 
HST observations was sufficient to examine their structure directly at optical wavelengths. An important 
finding was that the radial surface density profiles appear to be abruptly truncated at some outer radius, 
perhaps due to external effects from the surrounding Hll region and dense cluster (M096), and more 
detailed examination of this transition zone should lead to a greater understanding of the evolution of 
disks in harsh environments. 

The discovery images were obtained over a relatively narrow wavelength range (5007-6585A), and 
further insight should be possible through HST observations at shorter and longer wavelengths. In the 
blue/near-UV (~ 2000-4000A), the spatial resolution approaches ~15AU, while increased dust opacity 
at these wavelengths should also allow more tenuous structures to be traced to larger radii. Conversely, 
the considerable reduction in dust opacity at near-IR wavelengths should allow us to trace structures to 
smaller radii, albeit with commensurately poorer spatial resolution. Consequently, we are conducting 
follow-up HST studies from the near-UV to near-IR (0.3-2.5/im), and in the present paper, we report 
preliminary near-IR observations using NICMOS of one silhouette disk, Orion 114-426. The largest 
of the sample at ^1000 AU diameter, this disk is seen near edge-on, and while the central star is not 
directly visible at optical wavelengths, its presence is betrayed by two polar nebulosities believed to be 
illuminated by it. 
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2. 



Observations 



A comprehensive General Observer program (McCaughrean etai: GO 7367) studying the Orion 
silhouette disks with NICMOS, STIS, and WFPC2 is being carried out during HST Cycle 7. Early 
Release Observations using NICMOS were subsequently proposed by the Instrument Development Team 
(Erickson etai: SM2/ERO 7114) for scientific verification and media use. Due to this overlap, the ERO 
data were reduced and analysed collaboratively, resulting in studies of 114-426 (presented here) and of 
the 182-413/183-419 field (Chen etai. 1998). 

NICMOS observations of the 114-426 field were obtained on 19 April 1997 during the Servicing 
Mission Orbital Verification following installation in the HST. Images were taken through broad-band, 
narrow-band, and polarimetric filters between 1 and 2.1/xm as summarized in Table 1. Data reduction 
combined standard ground-based near-IR imaging techniques with parts of the NICMOS calibration 
pipeline. Multiple read-outs combined with multiple positions on the sky were used to reject cosmic-ray 
events; electronic offsets were removed with on-orbit dark images; quantum efficiency variations were 
removed with flat fields taken on-orbit where possible, otherwise from ground tests. Finally, mosaics were 
made registering the multiple images using stars or HST pointing information. Detailed photometric 
calibration was not attempted, but ground-based near-IR magnitudes for stars in the field were used to 
calibrate within ±0 m 2. 

Despite integration times significantly shorter than those planned for the GO program, important 
preliminary results were nevertheless obtained from the narrow-band imaging against the bright Paa 
background at 1.87/um, broad-band imaging at 1.1 and 1.6^m, and the polarization imaging at 2.0/im. 
The three polarizer position images were combined to form a 2/im continuum image, but due to 
remaining uncertainties in the correct analysis techniques for NICMOS polarimetry and incomplete 
on-orbit polarization calibration, the polarization results themselves are deferred to a future paper. 
The remaining narrow-band images did not provide useful additional information and are not further 
discussed. 



The highest S/N images of the silhouettes obtained by M096 were through a narrow-band Ha 
(A6565A) filter, admitting the full emission line flux from the bright Orion Nebula Hll region, while 
minimizing continuum emission from the central stars, or in the case of 114-426, its polar lobes. The 
brightest near-IR counterpart is the Paa line at 1.87/zm, which cannot be detected from the ground due 
to atmospheric absorption. For typical Hll region ionization parameters (10 4 K, 10 4 cm~ 3 , Case B) and 
Ay~l m foreground to the nebula, the detected photon flux at Paa should be ~60% of that at Ha: the 
brightest equivalent line available to ground-based observers (Br7 at 2.16/im) would be a further factor 
of ten fainter (Ostcrbrock 1989). 

The Paa 1.87/im image of 114-426 is shown in Figure 1 with the Ha (A6565A) image from M096. 
The S/N in the Pa image is poor (^5:1) since the integration time was short (288 sec), and the NIC1 
image scale of 0.0432 arcsec/pixel over-resolved the 0.19 arcsec FWHM diffraction-limited resolution 
of the telescope at 1.87 /im. Nevertheless, the silhouette is clearly detected, allowing a preliminary 
measurement of its size. The data were binned by a factor of two to better match the appropriate pixel 
size (i.e., 2 pixels per FWHM) and then averaged across the minor axis. The resulting ID major axis 
profile had high enough S/N to show the two ends of the disk as sharp dips separated by 1.8 arcsec. As 
discussed in detail by M096, the apparent size and structure of a silhouette disk is a convolution of its 
real form with the instrumental point spread function, and following M096, we adjusted the parameters 
of a model edge-on disk convolved with a model HST+NICMOS PSF calculated using the Tiny Tim 
software (Krist & Hook 1997) until the major axis length was reproduced. The resulting best-fit model 



3. 



Results 



3.1. Narrow-band imaging in Paa 



HIGH-RESOLUTION NEAR-INFRARED IMAGING OF ORION 114-426 



3 



disk has a major axis size of ~800AU at 1.87/xm, ^20% less than the 1012 AU measured at 0.6/.tm 
(M096). The same procedure used on the 2/zm continuum image (Section |3.2[ ) yielded the same result 
to within 5%. Finally, we verified the overall procedure by degrading the high S/N [OIII] image from 
M096 to the same spatial resolution and S/N as the Pact image, then performing the same fitting process, 
before retrieving the correct size for the disk at optical wavelengths. 



3.2. Continuum imaging of the polar lobes 

The optical continuum image of 114-426 showed faint polar lobes, interpreted as reflection nebulae 
of tenuous dust above and below the plane of the disk, illuminated by the otherwise unseen central star 
(M096). Similar reflection nebulae are seen above and below the plane of an edge-on disk in the HH30 
system (Burrows etal. 1996). The wavelength dependent morphology and polarization structure of the 
lobes in 114-426 should allow us to probe the underlying form of the disk, the geometry of polar cavities, 
the grain size, and scattering function. 

The near-IR broad-band continuum images, along with the F547M image from M096, are shown 
in Figure 2 in grayscale and contour forms. As the wavelength increases, three effects are seen. First, 
the initially fainter SE polar lobe increases in brightness until it equals then outshines the intensity of 
the NW lobe. The peak intensities in the lobes (after background subtraction) are in the ratios 7.3:1, 
2.2:1, 1.2:1, and 0.85:1 for the NW:SE lobes, at 0.57, 1.1, 1.6, and 2.0^m respectively. Second, the 
nebulae move closer together, as the reduced extinction allows us to probe closer to the disk midplanc: 
the separations of the peak pixels in the two lobes are at 0.64, 0.43, 0.40, and 0.32 arcsec at 0.57, 1.1, 
1.6, and 2.0/im respectively. Third, the nebulae appear to flatten from conical to slab-like. 

These features can be compared to model disks and envelopes {e.g., Lazareff, Pudritz, & Monin 
1990; Whitney & Hartmann 1992 [WH92], 1993 [WH93]; Fischer, Henning, & Yorke 1996 [FHY96]). 
The general broad fan shape and increasing flatness of the nebulae are best reproduced by model SH 
of FHY96, i.e., a thin, massless disk (in comparison to the central star) in Keplerian rotation, with an 
envelope and broad polar cavities. Models with thick, massive disks have more polar material and result 
in images with too much elongation perpendicular to the disk. The same is true of models with just a 
narrow, cylindrical polar hole rather than a broader, so-called "streamline" cavity (WH93). 

Since the central star in 114-426 is not seen, the disk must lie within a few degrees of edge-on, as the 
thin disk of the SH model does not occult the central star unless this condition is met (FHY96). This 
degree of alignment is also argued for on the grounds that the two lobes have nearly equal brightness in 
the near-IR. The asymmetry between the two lobes in the optical continuum (Figure 2a) can probably 
explained by asymmetries in the outer, more diffuse parts of the envelope, perhaps due to external effects 
in the H II region. 



4. Discussion 
4.1. Disk structure 

The disk appears ~ 20% smaller at 1.87/xm than at 0.6/zm, and thus we are clearly resolving structure 
in its outer parts. In order to understand the implications of these observations, we need to examine 
the theoretical expectations. For a thin, massless, Keplerian disk that is hydrostatically supported and 
vertically isothermal (Shakura & Sunyaev 1973; Pringle 1981; Lazareff etal. 1990), the density p as a 
function of radius, r, and height above the midplane, z, is: 



p(r, z) 



-15/8 



Pd 



exp 



4 \h{r), 



(1) 
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where pd is the midplane density at the outer disk radius, rd, and h(r) is the disk scale-height: 

9/8 

Jd 

where Zd is the scale-height at r d . The disk surface density is then: 

-3/4 

Jd ' 



h(r)=z d -) (2) 



E(r) = 2p d z d I — ) (3) 



For typical Orion Nebula dust grains (R=5; Cardelli, Clayton, & Mathis 1989), the extinction at 
1.87^m is one sixth of that at 0.6/mi, and thus achieving the same effective optical depth requires six 
times higher column density at the longer wavelength. For a face-on disk with unity optical depth at 

0. 6/iin at 506 AU radius, Eq.|| shows that the equivalent optical depth at 1.87^tm would occur at 46 AU, 

1. e., the disk would appear much smaller in the near-IR. For an edge-on disk, the calculation is harder: 
we are probing the midplane density of the disk, not the surface density, and the line-of-sight through 
the disk at a given "impact parameter" (i.e., the distance off-center) integrates over different densities at 
different radii. Assuming the disk is truncated at some outer radius, we can integrate the total column 
density through the midplane as a function of the impact parameter, a. We have calculated profiles for 
disks in which the midplane density scales as r a , with a = — 1, —2, and —3 (Figure 3). These values were 
chosen since simple analytical integrals exist in each case, but they also closely correspond to plausible 
disk models: a = —9/8 ~ —1 would yield a surface density independent of radius; a = —15/8 ~ —2 
yields the canonical S(r) oc r -3 / 4 ; a = —21/8 ~ —3 yields S(r) oc r -3 / 2 , a commonly assumed density 
law for more massive circumstellar disks (cf. Adams, Shu, & Lada 1988). 

Rough power laws can be fit to the inner section of each curve in Fig. 3, with the column density 
increasing as a 1 , and 7 ~ —1/4, —1, and —2 corresponding to a — —1, —2, and —3 respectively. In the 
canonical case of a = —2 ~ —15/8, a disk with unity optical depth at 506AU at 0.6/im would have 
the same optical depth at 84 AU at 1.87/mi. Thus, even though an edge-on disk will shrink less with 
wavelength than a face-on one, it is clear that for 114-426, the observed decrease in size of only 20% 
(i.e., — O.ldex) for an increase in density of six (i.e., 0.78 dex) would require 7 ~ —8, and thus a ~ —9 
and S(r) oc 7'~ 8 , quite inconsistent with conventional understanding of disk structure. 

The exception comes near the edge of a disk, where the integrated column density rises sharply as 
the path length increases most rapidly. Indeed, Fig. 3 shows that in the a = —2 case, the density does 
increase by the required factor of six within 80% of the outer radius, as observed. Thus, the observed 
shift between the optical edge at 506 AU and the infrared edge at 400 AU could still be consistent with 
the midplane density following the canonical power law for a thin, massless disk, as long as the disk is 
physically truncated somewhere near the optical edge. As discussed by M096, it seems plausible that 
the disk may be truncated due to physical processes present in the Orion Nebula Hll region and/or 
Trapezium Cluster, and it may be possible to ascertain which particular process is responsible from our 
future HST observations, which will trace the radial structure of 114-426 over a full order of magnitude 
in limiting column density. 



4.2. The disk mass revisited 

The mass of the 114-426 disk remains unknown: in particular, does it exceed the ~0.01M Q of the 
minimum mass solar nebula (MMSN), i.e., is it massive enough to form a planetary system similar to 
our own? O'Dell & Wen (1994) and M096 discuss how the mass of a silhouette disk can be roughly 
estimated, noting that the minimum intensity seen towards it is typically ~10% of that of the background 
Hll region. Assuming this to be background emission attenuated by the disk, a line-of-sight column 
density is calculated for each point in the silhouette, and summing over the whole disk area, a total mass 
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is estimated. However, virtually all of this "disk light" is an artifact, as the the background emission 
is blended into the disk by the instrumental PSF (M096). Thus the usefulness of the "transmission 
technique" is diminished, but it at least provides an absolute lower-limit mass, an important counterpoint 
to the results from millimeter interferometry which are, as yet, generally provide only upper-limit masses 
for disks in the Trapezium Cluster (Mundy, Looney, & Lada 1995; Lada etal. 1996). 

M096 used the transmission technique to estimate lower-limit masses for the six silhouette disks 
based on their optical images. The masses were small, typically several orders of magnitude less than 
the MMSN. The most massive was 114-426 at 0.002 M©, but we must report here that an error was 
made in calculating the area and thus mass of 114-426 in M096: its minimum mass should be revised 
downwards to 2 x 10~ 4 M©. (The masses for the other disks given in Table 2 of M096 remain unaffected 
by this error). 

In principle, the transmission technique could be applied to our new near-IR images of 114-426, 
resulting in a significantly higher estimate of the minimum mass. Since the disk is only marginally 
smaller at 1.87/xm, the reduced near-IR extinction should lead to a higher mass, on the order of the 
factor of six in Ay/Ai ^^. However, the present images are not good enough to make such estimates: 
the Paa image has very poor S/N, making a reliable measure of the disk/background flux ratio impossible, 
and while the continuum images have much higher S/N, the silhouette is contaminated by the reflection 
nebulae, making it difficult to measure the underlying disk/background flux ratio and the true silhouette 
area. Such calculations must wait until higher-quality Paa data are available. 

In the interim, alternative approaches to estimating the disk mass can be taken. The non-detection 
of the central star in the POL 2/im continuum image can be used to set a lower limit on the total column 
density through the disk midplane, depending on the intrinsic brightness of the star. The integrated 
flux from the two reflection nebulae is K—14™5, which clearly only represents a fraction of the flux we 
would see from the central star if it had no disk. Some estimate of this fraction can be made from 
WH92 and WH93, who tabulate F/F+ (the total flux observed as a fraction of the flux that would be 
seen if the star had no disk or envelope) for a variety of models. While F/F* ranges from 0.01-10% 
depending on the disk and envelope structure, inclination angle, grain properties, and total opacity, 
~1% is typical for models inferred for 114-426 (i.e., near edge-on flared disk, e.g., WH92 models 3,4 
with n < 0.3, perhaps with a tenuous envelope with broad polar cavities, e.g., WH93 models 5-8). The 
inferred intrinsic brightness of the central star is then ^T^9™5, corresponding to a 1 Myr old, ~1.5M© 
star at 450 pc (D'Antona & Mazzitelli 1994). Since the star is undetected in the 2/im image to a 5<r 
point source detection limit of iiT=18 m 8, the implied line-of-sight extinction towards it is >9™3 at 2pm, 
or Ay > 60™. Assuming the disk is exactly edge-on, the corresponding column density through the disk 
midplane to the central star is > 1.1 x 10 23 cm~ 2 (assuming lAy = 1.9 x 10 21 cm~ 2 ). 

The disk outer radius, rg = 506 AU, is known from the data of M096. The same data can be used 
to determine Zd, the scale-height at r^: we fit a ID average minor axis profile across the [OIII] image 
of 114-426 (Fig. 4 of M096) with a Gaussian characterized by Eq. |[ yielding Zd — 72 AU. Finally, we 
assume the disk inner radius, Vi = 1 AU: since the surface density increases towards the center with a 
power < 1, our results are not too sensitive to this assumption. Integrating Eq. [l] (with a = —15/8) 
through the midplane, we calculate the midplane density at the outer edge, pd > 5.8 x 10 4 cm -3 . Then 
integrating Eq.|| between and n, the total number of particles (HI + 2H 2 ) is > 3.4 x 10 52 , yielding a 
mass estimate of > 5.7 x 10 28 g or > 10 M 8 . 

Although it is clear that very little mass is required to render the central star invisible, the disk 
mass is likely to be significantly greater. Using the analytical solution for the total column density seen 
through the midplane of an edge-on disk in the a — —2 ~ —15/8 case (Section 4.1), and assuming 
Pd = 5.8 x 10 4 cm~ 3 at an outer radius of 506 AU, we obtain values ~ 2 x 10 20 cm~ 2 or Ay ~ m l just 
inside the edge. Yet the disk must have a significantly higher column density there to be seen as an 
edge: the results of M096 show that the attenuation at the edge of the disk is at least 85%, equivalent 



G 



McCAUGHREAN ET Ah. 



to Ay > 2 m , in turn implying pd > 1.2 x 10 6 cm~ 3 . This is a plausible value if (for example) the location 
of the disk edge is set by pressure balance with the surrounding H II region: assuming a temperature 
and density for the latter of 10 4 K and 10 4 cm -3 and an outer disk temperature of 10-100 K, densities 
at the disk edge of 10 6 -10 7 cm~ 3 would be predicted. 

Inserting pd > 1.2 x 10 6 cm~ 3 into Eq. [| along with our other derived disk parameters yields a disk 
mass of > 1.1 x 10 30 g or > 5 x 10~ 4 M Q . While still considerably less than the MMSN, this is a lower 
limit, and the disk might be more massive. However, our present simple calculations based on lower-limit 
extinction measurements alone cannot be pushed much further. Burrows etal. (1996) have shown for 
HH30 that more realistic mass estimates can be made by fitting the structure of the polar reflection 
nebulae, in particular the width of the dark lane of obscuration between them, and such an approach will 
be taken with 114-426 when our higher S/N and resolution Cycle 7 HST data are available. Ultimately, 
the most meaningful estimate of disk mass must come from measurements of optically thin tracers at 
millimeter wavelengths: such observations have been obtained for 114-426 and will be discussed in a 
future paper. 

5. Conclusions 

We have presented preliminary NICMOS images of the Orion 114-426 circumstellar disk. The disk 
appears to be roughly 20% smaller at 1.87^m than at at 0.6/um, which we attribute to a standard 
radial power-law density distribution inside a heavily truncated edge to the disk at ^500AU radius. 
The general wavelength-dependent morphology of the polar lobes is consistent with models of a thin, 
massless disk surrounded by a tenuous envelope with broad polar cavities, and their flux was used as 
an indirect way of estimating the mass of the central star to be 1.5 M®. The non-detection of the star 
at 2/im was then used to estimate the minimum total column density through the disk midplane, and 
thus a minimum disk mass of > 10 Mq, assuming that it follows the standard form for a massless disk 
in Keplerian rotation. The hard edge of the disk was used to improve the lower-limit estimate under 
similar assumptions, resulting in a mass > 5 x 10 _4 M Q . 
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TABLE 1. Observational log 



Camera 


Filter 


Integration 


Notes 


(a) 


(b) 


(c) 




NIC2 


POLOL 


384s 


1.9-2.1 /_im, degree polarizer 


NIC2 


POL120L 


384s 


1.9-2.1 /_im, 120 degree polarizer 


NIC2 


POL240L 


384s 


1.9-2.1 /^m, 240 degree polarizer 


NIC2 


F212N 


384s 


2.121 fim, 1% wide, H 2 v=l S(l) 








line 


NIC2 


F215N 


384s 


2.150 //m, 1% wide, continuum 


NIC1 


F110W 


384s 


0.8-1.4 fim 


NIC1 


F160W 


384s 


1.4-1.8 fim 


NIC1 


F164N 


288s 


1.644 fim, 1% wide, [Fell] line 


NIC1 


F166N 


288s 


1.650 fim, 1% wide, continuum 


NIC1 


F187N 


288s 


1.875 /j,m, 1% wide, Pa line 


NIC1 


F190N 


288s 


1.900 fim, 1% wide, continuum 



(a) NIC2 image scale 0.0757 arcsec/pixel, field-of-view 19.2x19.2 arcsec; 
NIC1 image scale 0.0432 arcsec/pixel, field-of-view 11.0x11.0 arcsec 

(b) Standard HST nomenclature for filter names used 

(c) Total integration time, in each case split into three sub-images at slightly 
dithered positions, each taken using NICMOS MULTIACCUM integration 
sequences for improved signal-to-noise and cosmic ray rejection 
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i i i I i i i i I i i i i 
(a) F656N: Ha 




(b) F187N: Paa 



Fig. 1. — Emission-line images of the 114-426 silhouette disk, (a) Ha 6565A, (b) Paa 1.87/xm. Both 
images have been magnified and rotated to a common scale and orientation. Each image is 4 x 4 arcsec 
in size or 1800 x 1800 AU at 450 pc. Major tick marks are at 1 arcsec intervals. The intensity scaling is 
arbitrary. 



10 



McCAUGHREAN ET Ah. 





(c) F160W: 1.4-1.8 /im (d) POL L: 1.9-2.1 /zm 



Fig. 2. — Continuum images of the 114-426 disk and polar lobes in grayscale and contour form, (a) 
F547M (0.52-0.58/xm), (b) F110W (0.8-1.4/xm), (c) F165W (1.4-1.8/zm), (d) intensity image derived 
from the POL0L/POL120L/POL240L polarization images (1.9-2.1/zm). All images have been magnified 
and rotated to a common scale and orientation. Each image is 2.86x 1.71 arcsec in size or ~ 1290x770 AU 
at 450 pc. Major tick marks are at 1 arcsec intervals. The absolute intensity levels are arbitrary, but the 
contours are spaced by 0.1 dex. Contours below the local mean background tracing the silhouette disk 
are drawn as dashed lines. Note that the fainter "side-lobes" seen above and below the main nebulae in 
the POLL image are artifacts resulting from the diffraction PSF of the HST+NICMOS. 
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Fig. 3. — Model midplane column density profiles for an edge-on disk with an outer edge at unity radius. 
The total integrated column density (in arbitrary units) is shown as a function of the impact parameter, 
i.e., the distance off-center at which the column is being measured. The three curves correspond to 
underlying radial power-law density distributions of the midplane density with p oc r a , with a = — 1, — 2, 
and —3. See text for further details. 



